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ABSTRACT 


This raport documents modeling end control system design for the 
Large Space Systems Tschnology (LSST) Reference Platform. The LSST 
Reference Platform consists of e central bus structure, solar panels, and 
platform arms on which a variety of experiments may be mounted. Simple 
structural models and classical frequency domain control system designs 
are developed. The report shows that operation of multiple Independently 
articulated payloads on a single platform presents major problems when 
subarc second pointing stability Is required. Experiment compatibility 
.will be an Important operational consideration for systema of this type. 
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I. INTRODUCTION 

This report documents modeling and control system design for the 
Large Space Systems Technology (LSST) Reference Platform, The LSST 
Reference Platform consists of a central bus structure, solar panels, 
and platform arms on which a variety of experiments may be mounted. 
Figure 1-1 shows one possible configuration. 

The objective of this report Is to Identify major control problems 
associated with precision pointing of experiment payloads. Simple 
structural models and classical frequency domain control system designs 
are used. 


A major result of this report Is to show that operation of multiple 
independently articulated payloads on a single platform will present 
major problems when high performance is required. Experiment compatibility 
will be an Important operational consideration for any systems of this 
type. 


The report which follows has six sections and six appendices. 

Section II defines the structural model for the platform. Mode shapes 
and frequencies are presented and controllability/observability issues 
discussed. Section III defines control system requirements. In Section 
IV the control system design is developed. Rate plus position feedback 
controllers are used. Control system design is evaluated in Section V, and 
in Section VI a summary is given along with comments on future study options. 

The Appendices document the model development for the platform 
structure. The model is developed in a building block fashion using 
finite element techniques. In Appendix A a model for the flexible 
piatform cross arms is given. In Appendix B the flexible solar panels 
are added. Appendix C presents transfer functions for the model of 
Appendix B. Two-hinged experiments are added to the model in Appendix D, 
and in Appendix E a technique for reshaping rigid body modes is given. 
Appendix F presents time response plots in support of Section V. 
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II. STRUCTURAL MODEL 

This seccion presents mode shapes for the LSST Reference Platform. 
Techniques developed in Appendix E are used to reshape the rigid body 
mode shapes so that they are easier to Interpret. Some comnents on 
controllability and observability are then made. 

A. CONFIGURATION 

Figure 2-1 shows a model of the spacecraft configuration. This 
configuration Is the same as that presented iu Appendix D. The variables 
1 ■ 1,2 ,,.., 6 are translational degrees of freedom; 6^ 1 1,2,3 and 

Yj^, Y3 are rotational degrees of freedom. Vj^, V2, V^, 6^, 62, and 63 are 
associated with the platform arms, V^, V3, and with the solar panels 
and Y^, Y3 with the experiment packages. The central bus mass and inertia 
are represented by M2, I2, and I5(M2 ■ « 1/2 total bus mass, 

I2 ■ I3 * 1/2 total bus Inertia). The bus Is assumed to be rigid. 

The experiment packages are hinged with respect to their bases. These 
bases have masses , M3 and inertias I^, and are connected to the 
central bus by massless beams of length L and flexural rigidity El . 

^ SL 

The experiment packages have masses , M^ and Inertias about their 
hinge point of 1^^, I3. The distances from the hinge point to the 
experiment package mass centers are L^, L3. The experiment packages 
and their connecting arms are assumed to be rigid. 

The solar panels are represented by M^, 1^, M^, (M^ “ ^6* ^4 " ^6^* 

The masses M4 and N^ are connected to the central bus by massless beams 
of length and flexural rigidity £1^. 

B. PARAMETER VALUES AND NATURAL FREQUENCIES 

Table 2-1 lists two sets of parameter values used for simulation 
purposes. Parameters for the bus and solar panels are similar to those 
of Appendix B. Parameters for the experiment packages approximate those 
of SIRTF* . As can be seen frexn Table 2-1, the two data sets are Identical 
except for the length of the platform arms. 

^Shuttle Infrared Telescope Facility (See Aviation Week, Sept. 15, 1980). 
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Table 2-2 Hats the natural frequencies for the 6 elastic modes of the 
model. The primary difference between set 1 Md set 2 is the change in 
frequency of the platform arm modes which results from shortening the 
arms per Table 2-1. 

C. MODE SHAPES 

The mode shapes for the S rigid body modes and first 4 elastic modes 
are sketched in Fig. 2-2 and Fig. 2-3. These mode shapes do not change 
qualitatively for the different data sets. 

Numerical values for the first 9 mode shapes for the parameters of 
Table 2-1, set 1, are given In Table 2-3 for the coordinates of interest 

^ 2 * ^ 3 * ^ 1 * ^ 3 ^* coordinates are associated with activation 

and/or sensing. The center of mass for the solar panels (M^, M^) and 

platform arms (M^, M 2 , M^, M^, M^) has been appended to each mode shape. 
Note that these mass centers are not associated with any physical point 
of the structure. 

The rigid body mode shapes of Fig. 2-2 and Table 2-3 are more 
complicated than they need to be. For this reason, the technique developed 
in Appendix E was applied to obtain more pleasing sh^es. The resulting 
shapes are presented in Fig. 2-4 and Table 2-5. 

In order tc obtain the rigid body shapes the following steps were 
used. First it was desired to eliminate translation of the solar panel 
center of mass from all but one mode. The computer algorithm selected 
the mode shape with the largest value of solar panel center of mass 
translation (mode 3 of Table 2-3) and using this mode shape as described 
in Appendix E it eliminated this coordinate from modes 1, 2, 4, and 3. 

This resulted in five new mode shapes. 

Next the algorithm was told to eliminate the arm center of mass 
translation from all but one of ihe modes 1, 2, 4, or 5. Mode 1 was 
selected and four new mode shapes found. 

Next, the algorithm was told to eliminate the bus rotation angle ( 62 ) 
from all but one of the modes 2, 4, or 3. Mode 3 was selected and three 
new mode shapes were found. 
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Finally, tha algorithm waa told to allmlnata tha axperiment rotation 
angle (y^) from all but one of the modea 2 or 4. Mode 4 waa aelected and 
two new mode ahapea reaulted. 

It should be noted that if a different set of coordinates or a different 
ordering had been selected, then different rigid body mode shapes may 
have resulted. The reason for selecting the mess center coordinates will 
be more clear when controllability/obasrvability issues are discussed 
later in this section. 

D. CONTROLLABILITY 

The controllability of various modes can be seen as follows: 

The equations of motion for the platform are 

Mx + Dx + Kx ■ Bu 

where 

M, D and K - 

mass, damping, and stiffness matrices for the structure, respectively. 

X ■ physical coordinate vector (dimension n) 
u ■■ control vector (dimension m) 

B * control distribution matrix (n x m) 

In modal coordinates this equation becomes 

•• T • . 2 T 

D^q + A q*^ Bu 

where 4 (the modal matrix) satisfies 

x ■ ♦ q 

M ♦ - 1 

n 

K ♦ • - diag. (A^, A^, ..., A^) 

A^ ■ i^^ natural frequency 

T 

If we assume modal damping, then « D4 - D is diagonal. Hence, we see 

T 

that the matrix product 4 B determines the ability of the control u to 
influence the modal dynamics. 
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The control forces (torques) applied to the platform are given by 
F - Bu 

where 

u - 

<■ actuator torque vector (T2 Is central bus torque, and 
are experiment torques) 


B 


0 

0 

0 

0 

0 

2.07E6 

0 

0 

0 

0 

0 


0 

0 

0 

-4.8E4 

0 

0 

0 

0 

0 

4.8E4 

0 


0 

0 

0 

0 

-4.8E4 

0 

0 

0 

0 

0 

4.8E4 


o Control distribution matrix 

Notice that B Includes torque scaling for each actuator. Notice also that 
the actuators which articulate the experiment systems (T^ and T^) apply 
a torque to the experiments and an equal and opposite toique to the platform 
arms at the point of attachment. 

T 

Table 2-6 gives the matrix product $ B for the eigenvectors of Tables 2-4 

T 

and 2-5, Here, only the first 9 modes of 4 are Included, so that 4 B Is 
a 9 X 3 matrix. 


We can determine the ability of an actuator to Influence each mode 

T 

directly from the elements of 4 B. For example, mode 1 (arm translation), 
mode 3 (solar panel translation) , and mode 6 (solar panel symmetric bending) 
cannot be Influenced by any of the three torque actuators and hence these 
modes are uncontrollable. Modes 4 and 8 can be Influenced by T^ and T^ 
but not T 2 . 
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It should be noted that for distinct eigenvalues, the ability of a 
force to Influence a mode Implies controllability of that mode (see 
Ref . 2>1) . For modes associated with repeated eigenvalues , however , this 
Is not the case. Since there are five rigid body modes, then nt least 
five actus tors are required to provide controllability of these modes. 
Clearly torque actuators cannot move the spacecraft mass center, so any 
mode shape having a spacecraft mass translation Is not controllable. In 
Table 2-3 all of the rigid body modes have mass center translation. In 
Table 2-5 only modes 1 and 3 have CM translation, so that modes 2, 4, 
and 5 form a controllable subset of rigid body modes and modes 1 and 3 
form an uncontrollable subset. 


E. OBSERVABILITY 

The sensors we will consider measure positions and rates of 62 > Yj^> 
and Since we have 


X “ 4 q 


the outputs are y = Cx where 


X . [Vj Vj V3 6^ 63 63 V4 Vj Vj n Y3I' 


c 


Now 

y 


000 

000 

000 

Cx = C$q 


0 0 1 
0 0 0 
0 0 0 


0 0 0 0 0 
0 0 0 1 0 
0 0 0 0 1 


so determines Che Influence of a given mode on the sensor outputs. 
Table 2-7 gives Ct for the eigenvectors of Tables 2-4 and 2-5. 


From Table 2-7 we see that modes 1, 3, and 6 do not influence the 
system outputs, and hence are unobservable. A comparison of Table 2-7 
with Table 2-6 shows that actuator and sensor Influence results are Che 
same (e.g. if is influenced by Tj^, then it influences the sensor output 
in Y^. Similarly for T 2 and 62 , and for T^ and Yj*)* 
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For distinct eigenvalues, the ability of a mode to Influence a sensor 
output Implies observability of that mode. For repeated eigenvalues, 
however, this is not the case. A minimum of five sensors would be required 
for observability of the rigid body modes of Table 2-3 or Table 2-5. 

The three rotation sensors have no capability to sense translation of 
either the platform arms or the solar panels. In Table 2-3 modes 2, 4, 
and 5 form an observable subset of rigid body modes, and modes 1 and 3 form 
an unobservable subset. 

Since modes 1, 3, and 6 are neither influenced by the actuators, nor 
do they Influence the sensor outputs we may eliminate them frtnn further 
consideration in control system design. 


Table 2-1. 

Parameter 
Platform Arms 

El (N-m^ X 10^) 

& 

(k« X loh 

M 2 (kg X 10^) 

M 3 (kg X 10^) 

(kg-m^ X 10^) 

1 2 (kg-m^ X 10^) 

1 3 (kg-m^ X 10 ^) 

Solar Panels 
Ljj (m) 

EIjj (N-m^ X 10^) 

M^ (kg X 10 ^) 

M 3 (kg X 10 ^) 

M^ (kg X 10^) 

I 3 (kg-ro^ X 10 ^) 

Experiment Packages 
L^l (m) 

Le3 (m) 

Mgj^ (kg X 10 ^) 

M ^3 (kg X 10^) 

(kg-m^ X 10^) 

I - (kg-m^ X 10^) 


Parameter Values 


Set 1 

Set 

15 

6 

20 

20 

0.3 

same 

6.0 

II 

0.3 

II 

0.1 

M 

50.0 

II 

0.1 

tl 


20.0 

same 

0.1 

It 

0.5 

II 

6.0 

It 

0.5 

II 

Included 
In I 2 

II 

3.6 

same 

3.6 

II 

3.3 

II 

3.3 

II 

48.0 

II 


A8.0 












Table 2-3. Rigid Body Mode Shapes Before Reshaping 


RIGID BODY MODES 


2 



CM (ARMS) 


6.14 E-3 


-3.46 E-3 

2.11 E-7 

-3.08 E-3 

3.82 E-7 

1.06 E-6 

j 

1.20 E-2 

3.30 E-4 

1.74 E-7 


7.23 E-6 

1.15 

E-3 

7.23 E-6 

1.15 

E-3 

7.23 E-6 

1.15 

E-3 


-4.59 E-3 
-3.85 E-3 


-3.22 E-13 
-8.38 E-4 



Table 2-4. Elastic Mode Shapes 


ELASTIC MODES 
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3.28 E-18 1.25 E-3 -3.04 E-3 j 3.12 E-3 


2.97 E-18 , 1.25 E-3 


3.04 E-3 3.12 E-3 


3.03 E-13 1.24 E-3 I 3.01 E-11 -4.14 E-3 


6.16 E-3 -2.57 E-3 

-6.16 E-3 -2.57 E-3 


•8.56 E-19 2.29 E-20 

2.13 E-11 -3.53 E-11 
























Coordinate 


Table 2-5. Rigid Body Mode Shapea After Reshaping 


RIGID BODY MODES 


3.28 E-10 

1 -5.28 

E-12 : 

1.66 

E-15 

1.01 

E-11 

1.16 

E-3 

3.68 E-10 

i 1.85 

E-11 

1.67 

E-15 

-2.25 

E-11 

1.16 

E-3 

3.38 E-10 

1 3.24 

E-14 i 

1.67 

E-15 

1.01 

E-15 

1.16 

E-3 


Yl 

-1.21 E-9 

— 
-2.10 E-11 

-5.28 E-14 

-5.40 E-3 

^3 

-1.34 E-9 

-5.18 E-3 

4.53 E-14 

1.55 E-3 

CM (S.P.) 

-2.38 E-15 

2.37 E-15 

1.20 E-2 

-1.24 E-15 

CM (ARMS) 

8.70 E-3 

-1.52 E-11 

-1.88 E-13 

2.02 E-11 


Table 2-6. ♦ B for the Eigenvectors of Tables 2-4 and 2-5 


T, 


1 

0 

0 

2 

0 

0 

3 

0 

0 

4 

0 

-2.59 E2 

5 

2.39 E3 

-2.62 E2 

6 

0 

0 

7 

2.56 E3 

-2.81 E2 

8 

0 

4.42 E2 

9 

-8.56 E3 

-2.73 E2 


0 


■2.48 

E2 

0 

7.43 

El 

2.63 

E2 

0 

•2.81 

E2 

■4.42 

E2 

■2.73 

E2 


Table 2-7. C4 for Eigenvectors of Tables 2-4 and 2-5 


f!nnr— j 




MODE 




1 1 

dlnate [ 1 

2 


4 

' 5 

Lll 

7 

8 ! 

9 


0 1.16E-3 0 1.24E-3 0 j-4.14E-3{ 

j 

5.40E-3|-4.29E-3 U -4.61E-3 6.16E-3 -2.57E-3| 

1.55E-3i-4.29E-3 0 -4.61E-3 -6.16E-3 -2.57E-3 































MODE 1: ARM TRANSLATION 



MODE 2: EXPERIMENT ROTATION, ANTISYMMETRIC 




MODE 3: SOUR PANEL TRANSLATION 


MODE 4: EXPERIMENT ROTATION, SYMMETRIC 


MODE 5: BUS ROTATION 





Fig. 2-2. Rigid Body Modes Before Reshaping 








Fig. 23. 


First Four Elastic Modes 


MODE 1: ARM TRANSLATION 


MODE 2: EXPERIMENT ROTATION 


MODE 3: SOLAR PANEL TRANSLATION 


MODE 4: EXPERIMENT ROTATION 


MODE 5: BUS ROTATION 




Fig. 2-4. 


Rigid Body Modes After Reshaping 


III. CONTROL SYSTEM REQUIREMENTS 


This section gives a brief statement of the controls, objectlvea, and 
requirements for the platform. 

A. OVERVIEW 

There are many requirements which attitude control systems for the 
LSST Platform have to satisfy. These Include a capability for: 

(1) stabilization and Initial acquisition of celestial references 

(2) maintaining a prescribed attitude In the presence of various 
disturbance Inputs 

(3) experiment 'olnting control 

(4) reorientinc the gross pointing direction to obtain a new viewing 
configuration or for reboost 

(5) reboost (or station keeping) to maintain a specified orbit 

Of these five requirements, only the second and third will be addressed in 
what follows. 

B. STABILITY REQUIREMENTS 

Requirements for attitude control Include both accuracy and stability. 
In this report we shall limit ourselves to a consideration of stability. 

In addition to pointing stability we will also consider linear acceleration 
levels, since some experiments (e.g. material processing experiments) are 
concerned with these. 

Table 3-1 lists requirements for pointing stability and acceleration 
levels for some of the more stringent missions being considered for the 
platform. Pointing stabilities In the range of .01 to 1 arc sec will 
almost certainly require sophisticated pointing systems. The shuttle 
Infrared telescope facility has pointing requirements In this range 
( .1 arc sec) . 

Pointing stabilities In the range of 10 to 1000 arc sec will require 
less sophisticated pointing systems. In this range many users would prefer 
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the basic platform control system to meet their stability requirements. 
The cryogenically cooled telescope has pointing requirements in this 
range (30 arc sec). 

Pointing stability requirements below .01 arc sec will probably 
require image motion compensation. This problem will not be addressed 
in this report. There are a number of users with pointing requirenwnts 
between 1 and 10 arc seconds. Most probably they will require pointing 
systems although there are some users that would like the basic platform 
to provide this pointing capability. The high energy gamma telescope has 
pointing requirements in this range (6 arc sec) . 


Table 3*1. Platform Requirements 


Requirement 

Specification 

pointing system stability 
platform stability 
platform acceleration 

.03 to 5 wrad (.01 to 1 arc sec) 
.03 to 3 mrad (10 to 1000 arc sec) 
.1 to 10 mm/s^ (10 ^ to 10 ^ g) 
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IV. CONTROL SYSTEM DESIGN 

Thi« ••ccion docuaents control syatoa design for the LSST Reference 
Pletform of section II. The purpose of this design effort was to obtain 
some simple control system designs which would be suitable for Identifying 
the major control problems associated with the LSST reference platform. 

No attenq>t has been made here to obtain designs having optimal performance 
In any sense. Classical frequency design techniques are used based on 
simplified rigid body models of the platform. Rate plus position feedback 
controllers are used with ideal actuators and sensors assumed. The bus and 
experiment package controllers are designed Independently. Before beginning 
the control system design, several control concepts are discussed. 

A. CONTROL CONCEPTS 

Three control concepts for the LSST platform will be discussed in 
what follows. The first technique is the simplest conceptually and provides 
Independent control for the bus structure and experiment packages. The 
second scheme allows for limited one-way infonuation exchange between the 
bus and experiment package control systems. The final concept provides 
limited, or unlimited two-way Information exchange between control systems. 

1. Independent Control 

Figure 4-1 illustrates the various control schemes symbolically. 
Only two experiment packages arc Illustrated, but in actuality there 
could be many more. Independent control Ls the simplest conceptually. 

**ch control system operates independently using its own set 
of sensors and actuators. 

2 . One-Way Interface 

With the one-way interface scheme the central bus can send 
Information to each experiment package. In this case the bus controller 
operates independently, but the experiment package controllers take 
into account the data from the central bus structure in some way so 
as to improve pointing accuracy. 
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The oneway Interface scheme could also be implemented in the 
other direction. That is, each experiment control system could send 
the bus controller information which would help it to minimize the 
effects of disturbance inputs Introduced by that experiment package. 

This might be as simple as a status bit indicating whether or not the 
experiment control actuators are torquing, or, it could be a more complete 
Information exchange. With this scheme the experiment package control 
systems would operate independently, but the bus controller would not 
be totally Independent. This one-way interface approach would make 
the bus controller design dependent on varying experiment package 
designs and as such may not be very practical. 

i. Two-Way Interface 

This method gives each control system access to Information 
available to the other control systems. This might be a partial or 
total Information exchange. With this scheme the controllers are all 
Interdependent to some extent. This dependence could be complete, so 
that In effect there is only one central controller, or it could be 
less complete allowing more independence of action. As t te number 
of experiment packages increases, this two-way interface could become 
rather complicated. The two-way interface also makes the bus controller 
design dependent on the various experiment package designs and as such 
may not be very practical. 

4. Sensors and Actuators 

As a minimum we assume that each experiment package has some 
means of torquing to maintain Its desired orientation. A two-axis 
(e.g. azimuth, elevation angle) torquing capability might be typical. 

In addition, each package may or may not have its own (angular) 
position and rate sensors. If an experiment package has position 
sensors , it may or may not have the ability to directly measure target 
position. 

The bus controller will have a torquing capability plus position 
and rate sensors. As a minimum, sensors and actuators will be located 


on the central bus structure, but additional actuators and sensors 
could be located on the flexible portion of the structure If 
necessary to obtain the desired performance. 

The Information exchange, then, can be one In which sensor and 
actuator Information Is exchanged, and/or one in which estimated 
state vector Information Is exchanged. 

5 . Design Approach 

The approach to designing a control system for the LSST platform 
should be in the order of conceptual complexity. First the Independent 
controls approach, then the one-way interface. The two-way Interface 
may or may not be a viable alternative and should be attempted only 
as a last resort. In fact any scheme which makes the bus controller 
dependent on the various experiment packages would be an operational 
nightmare and should be avoided If at all possible. 

The bus controller should be attempted first with all actuators 
and sensors on the rigid central structure. However, If performance 
improvement is needed, the location of sensors on portions of the 
flexible structure might be considered. Location of actuators on the 
flexible structure is much more difficult and should be attempted only 
if other approaches fail to provide the required performance. 

The remainder of this report considers only the simplest design approach, 
that being independent control systems for the central bus and each experiment. 
It will be af.sumed that each control system has its own torquers and its 
own position and rate sensors. All sensors will measure inertial position 
and rate. 


B. BUS CONTROLLER DESIGN 

A block diagram for the bus controller is shown in Fig, 4-2. Symbols 
in Fig, 4-2 are defined as follows 


6 ■ Commanded bus angle (nominally zero) 
c 

6 <■ Actual bus angle (relative to some Inertial coordinate system) 

T ■> Torque applied to bus by controller 

T. ■ Disturbance Torque applied to bus 
d 

J " Platform inertia 

A 

J “ Torquer gain (an estimate of J) 

K •* Position gain 
P 

■ Rate gain 

In Fi(j. 4-3 the block diagram has been redrawn in a format convenient for 

controller design purposes. Here is the rate to position gain also 

called K in what follows, 
rp 

The forward loop transfer function from Fig. 4-3 is 

A 

K J 

G(s) » 

JS 

The feedback loop transfer function is 

H(s) = K S + 1 

rp 

From this the closed loop transfer function is 


J. G£s} 

0 l+G(s)H(s) 

c 


A 

K J 

E 

A A 

JS +K JK S-HC J 
p rp p 



S +(K K 
P rp 


dk 

J/J 

'a a • 

J/J)s+K J/J 
P 


or 
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Now If J « J then we have 


fu) 

c 


S Nic K S4K 
p rp p 


The characteristic equation for the system is then 

S^+KK S+K - 0 
p rp p 

Next define the damping and natural frequency as q and respectively. 


so that 


» K 
n p 


2?u .= K |C 
n r rp 


From this we can determine and given desired values for and 5 


rp 


2 

(A) 

n 

(4-1) 

2C(o 

(4-2) 

n 

2C/«„ 

(4-3) 


Or, given Kp and K^p we can determine t and 


(*) » 
n P 


4 % 


ij K /*^r 
r p 


*5 ‘“n 
n rp 


For the purposes of this design, we will choose 


u s 0*01 Hz 
n 


0.0628 rad/s 


0.707 


2 ^ 


then 

K ■ ■ 3.94E-3 rad/s^ 

P n 

K ■ 88.6E-3 rad/s - . ^ 

r 

K » 22.5 s 
rp 

To complete the controller design we must have a value for J . Figure 2-1 
shows a sketch of the platform model. It Is not Immediately obvious how 
to calculate J, because of the hinged experiment packages. Suppose, 
however, that we assume that the experiment pointing angles (Yj^ and Y 3 ) 
remain zero. This Is equivalent to saying that the commanded (inertial) 
pointing angles are zero for Yj^ and Y 3 , and that the experiment controllers 
are functioning ideally. Then so far as the bus controller is concerned 
the mass of the experiment M ^ can be lumped with the platform mass , 
and similarly with and M^. Then J can be calculated as follows 

J = + M 3 + M^ 3 )(L^)^ + (M^ + M^)(L^)^ + ^2 ^5 

Then for the parameters of Table 2-1, set 1, we have 
J = (0.3 + 3.3 + 0.3 + 3.3)E3 (15)^ 

+ (0.5 4- 0.5)E3 (20)^ + 50E3 
= 2.07E6 kg-m“ 

For Table 2-1, set 2 data we have 

J = 7.2E3 ( 6 )^ + 1.0E3 (20)^ + 50E3 
= 0.709E6 kg-m*’ 

This completes the bus controller design. 

C. EXPERIMENT CONTROLLER DESIGN 

A block diagram for the experiment controllers is shown in Fig, 4-4. 
This figure is identical to Fig. 4-2 except for the additional commanded 

*The inertia of the experiment about its hinge point need not be included 
since we assume that the experiment does not rotate relative to an inertial 
reference frame. 


Inputs T and y • T and v are the connnanded torque and rate 
c c c c 

respectively . 

For the experiment controllers we will choose two values for 
the natural frequency: 

u = 1 Hz = 6.28 rad/s 
n 

and 

u = 0.1 Hz ■ 0.628 rad/s 
n 

In each case we will select 5 = 0.707. 

Now using Eqs. A-1 through 4-3 we can compute the gains K^, K^, 

and K . For a 1-Hz controller we have; 
rp 

K * 0 )^ ■ 39.4 rad/ s^ 

P n 

K = 240 ) * 8,88 rad/s 
r n 

K “ 24 / 0 ) “ 0*225 s 
rp n 

For the 0.1-Hz controller we have; 

K = 0 )^ « 0.394 rad/s^ 

P n 

K = 240 ) “ 0.888 rad/s 
r n 

K - 24/o)_ » 2.25 s 
rp n 


The value for J can be read directly from Table 2-1. 

2 

J = I I - » 48E3 kg-m 
el 

This completes the design of the experiment controllers. 


D. COMBINED DESIGN 

A block diagram for the combined system is shown in Fig. 4-5. The 
commanded inputs are shown 6^2 tke bus controller, for the experiment 
3 controller , and , T^ 2 . experiment 1 controller . For the 

purposes of design evaluation, only the experiment 1 controller has a non- 


Yci is not shown on this diagram, because for the simulations conducted, Vd 
was set to zero. A better design could have been achieved using a nonzero 
value for Y^ii this was not realized until after all simulations had been 
completed. 
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zero comnanded Input. The success of the control system design is to be 
measured by how closely the bus and experiment 3 controllers maintain 
6 2 and near zero In the presence of coomtanded slews of the experiment 1 

controller . 

E. CLOSED LOOP SYSTEM EIGENVALUES 

The open loop eigenvalues for the LSST reference platform are given 
In Table 2-2. Table 4-1 shows the system gains used for simulation purposes. 
Figures 4-6 and 4-7 show Che associated closed loop eigenvalues for these 
gain values, (uncontrollable/unobservable modes are not plotted unless 
noted otherwise) . 

Figure 4-6, Run F12, shows the closed loop eigenvalues for the 15-m 
platform arms when open loop torqulng is used for experiment 1. Note the 
2 closed loop poles at the origin. This is an indication that the rigid 
body mode associated with experiment 1 becomes uncontrollable when the 
feedback gains are set to zero. The bus controller frequency and damping 
are somewhat higher than the design values of .0628 rad/s and .707 damping. 

Also, the experiment 3 controller frequencies are somewhat higher than the 
design values, and critically damped. The elastic mode vibration 
frequencies all have low values of damping. 

Figure 4-6, Run FIO, shows the closed loop eigenvalues for 15-meter 
platform arms and 1.0-Hz experiment controllers. Note that the bus 
damping and frequency is close to the desired values of 0.707 and 6.28E-2 rad/s. 
The experiment control frequencies, however, are critically damped with 
frequencies of 5.0 and 59.0 rad/s. The elastic mode vibration frequencies 
(for the controllable modes) all have small values of damping. 

Figure 4-6, Run F13, is similar to Run FIO except that the eigenvalues 
associated with the experiment 3 controller are shifted lower because of the 
lower gains used. 

Runs F17 - F19 use the same estimate for the pltiUfuriu Inertia 
2 

(J *■ 2.07E6 kg-m ) as for Runs F12, FIO, and F13. The design value of 
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2 

0.709 kg-o was not used to demonstrate the effect of parameter error on 
system performance. The results are shown In Fig. 4-7. Note that in 
all cases closed loop system response remains stable. In fact, damping 
for the elastic modes of the platform Is actually Improved over many cases 
with 15-m platform arms. Based on the results shown In Figs. 4-6 and 4-7, 
we might say that the design Is adequate for the purposes of Identifying 
major control problems, although it is in no way an optimum design. 


Table 4-1. Control System Gains for Performance Analysis 


RUN 

ARM 

LENGTH 

CONTROL 

BANDWIDTH 

EXP. 1 GAINS 

EXP. 3 GAIN 


(m) 

w , 
nl 

(Hz) 

”n2 

(Hz) 

K 

P 2 

(rad/s ) 

(rad/s) 

K 

P 2 

(rad/s ) 

K 

r 

(rad/s) 

F12 

15 

* 

1.0 

0.0 

0.0 

39.4 

8.88 

FIO 

15 

1.0 

1.0 

39.4 

8.88 

39.4 

8.88 

F13 

15 

1.0 

0.1 

39.4 

8.88 

.394 

.888 

F17 

6 

* 

1.0 

0.0 

0.0 

39.4 

8.88 

F18 

6 

1.0 

1.0 

39.4 

8.88 

39.4 

8.88 

F19 

6 

1.0 

0.1 

39.4 

8.88 

.394 

.888 


ALL RUNS USE THE FOLLOWING DATA: 

* 2 

Bus Inertia est, J • 2.07E6 kg-m 

Exp. 1 and Exp. 2 Inertial est. J ■ 48.0E6 kg-m^ 

Bus Gains K ■ 3.94E-3 rad/s^ K ■ 88.6E-3 rad/s 
P ^ 


^ Open Loop Torqulng 



INDEPENDENT CONTROL: 


EXPERIMENT 

PACKAGE 

CONTROL 

SYSTEM 


BUS 

CONTROL 

SYSTEM 


EXPERIMENT 

PACKAGE 

CONTROL 

SYSTEM 


ONE-WAY INTERFACE: 



TWO^AY INTERFACE; 



Fig. A-1. 


Control Concepts for the LSST Reference Platform 










Alternate Bus Controller Block Diagram 











Fig, 4-5. 



Systea Configuration 
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V. CONTROL SYSTEM EVALUATION 

In this section the transient response performance of the combined 
system Is examined. Comments are also made on additional factors which 
can affect performance. 

A. TRANSIENT RESPONSE 

Computer simulations were performed to determine system performance. 
Table 5-1 summarises the results (for a definition of the symbols used 
In Table 5-1, see Figs. 2-1 and 4-5). All runs were made with a bus 
control frequency 0.01 Hz. The commanded torque was a 

square wave Input of -f 20 N - m, followed by an equal duration -20 N - m 
Input. The constanded torque durations were chosen to give the desired 
final commanded angle (y. ) at the end of the torqulng sequence assuming 
that the base point (M^) was Inertlally fixed. The commanded angle as 
a time function was taken as the second integral of torque according to 
Newton's law (again with the assumption that the base point was fixed). 

The experiment inertias of Table 2-1 were used (T^^^ ■ ^cl^* 

Table 5-1 presents the peak transient responses for the experiment 3 

pointing angle (y^) , the bus pointing angle ( 6 ^) and the bus acceleration 

(V 2 ) . A comparison of Tables 5-1 and 3-1 shows that y^ exceeds a stability 

requirement of 5 wrad by up to 70 times, and is 60 to 7000 times a 0.05-urad 

2 

requirement. Acceleration levels are within the 10-mm/s requirement but 
exceed a 0.1-mm/s^ requirement by 6 to 30 times. 62 performance is between 

2 and 5 mrad (400 to 1000 arc sec) . This is within the 5-mrad requirement 
but 40 to 100 times the 0.05-mrad requirement. From this we see that the 
pointing requirements of Table 3-1 cannot be easily met using controllers 
of the type considered here. 

A comparison of the six representative simulation runs illustrates 
several clear trends. It should be kept in mind that experiment 1 i^ 
creating the disturbance, and experiment 2 (fixed to the bus) and experiment 

3 are feeling the influence of the disturbance. Using a closed loop 

* 

In the simulations, the base point (Ml) was not Inertlally fixed, so that 
the open loop torque conmand was in error. 
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controller for experiment 1 (compare Rune PIO and F12) results In a smoother 

disturbance torque profile, and as such Improves pointing stability for 

•• 

both Y 3 and It also reduces acceleration levels at the bus (V 2 ). 

Increasing experiment 3 bandwidth (ccmpare Runs FIO and F13) Increases Y 3 
stability, but has an adverse affect on V 2 . Decreasing the arm length 
(compare Runs FIO and F18) Increases pointing stability of both Y 3 and 62 
and In most cases decreases acceleration level (V 2 ) . The reader Is referred 
to Appendix F for time plots of an example simulation (Run FIO) . 

B. ADDITIONAL FACTORS WHICH CAN AFFECT PERFORMANCE 

There are many factors not considered In this study which could have 
a major Influence on the absolute performance of the LS5T reference platform. 
Imperfect sensors and actuators, glmbal friction and flexibility, and more 
complex structural dynamics could all result In poorer performance than 
that presented here. On the other hand improved controller gain selection 
or more sophisticated controllers could improve the performance results. 

For example, base motion compensation could be added using an additional 
sensor for eacJi experiment package (an accelerometer) , and image motion 
compensation could be implemented for some types of experiment packages 
by the addition of another actuator (a secondary mirror drive for example) . 
For these reasons, the performance results of Table 5-1 should not be 
taken as absolutes in any sense. Yet Table 5-1 does indicate the difficulty 
which is faced when attempting to design a control system for a platform 
of this type. 


Table 5-1. Performance Results 




CONTROL BANDWIDTH 

PEAK RESPONSE 


ARM 

W 

W 



V, 


liNGTH 

"1 

"3 


»2 

2 

(mm/s^) 

RUN 

(m) 

(Hz) 

(Hz) 

(prad) 

(mrad) 

F12 

15.0 

* 

LO 

56.00 

5.11 

3.000 

FIO 

15.0 

LO 

LO 

7.59 

3.27 

a 964 

F13 

15.0 

1.0 

ai 

357.00 

3.25 

0608 

F17 

6.0 

* 

LO 

12.20 

3.10 

0848 

F18 

6.0 

LO 

LO 

3.13 

2.28 

0819 

F19 

6.0 

1.0 

ai 

254.00 

r ?6 

0 819 

*0PEN LOOP TORQUING; 





Y^ = 88 mrad, T - 20 N-m 
C] 

. w„ • 0.01 
"2 

Hz for all 

runs. 



r 


> i 



I 

I 

i 

i 



VI. SUMMARY AND FUTURE STUDY 


Several control problems for the LSST Reference Platform have been 
identified and quantified in this report. Perhaps the most important of 
these is that operation of multiple independent control systems on a 
single platform presents a major problem when high performance is 
required. Experiment compatibility will be an important operational 
consideration. Control system design is complicated by large shifts in 
structural parameters which occur as a result of variations in the number 
and location of experiments mounted on the platform. Structural vibration 
frequencies in the controller bandwidth further complicate the design 
problem. It has been found that convenMor l controllers miss performance 
requirements by a wide margin when these factors are taken into account. 

It should be noted that the lighter the platform is. Independent of 
its stiffness, the greater is the controller interaction problem. The 
best structural design solution, from this standpoint, would be to place 
the Queen Mary in orbit. This might post other problems, however, from 
socio-political -economic viewpoints . 

It should also be noted that problems with platform flexibility cannot 
be solved by simply making the platform arms more rigid. As Appendix F 
illustrates, elastic vibration of the solar panels is the dominant pointing 
error once the initial slew transient (of experiment 1) has subsided. 

Two control approaches have been identified for future study. The 
first approach is to add additional sensors and/or actuators to individual 
experiment controllers. Base motion and image motion compensation fall 
in this category. The second approach is to allow information exchange 
between controllers, particularly one-way exchange from the bus controller 
to the experiment controllers. The challenge is to develop controllers 
which can significantly reduce the controller interaction problem and at 
the same time reduce controller sensitivity to structural parameter 
variations. 


PRECEDING PAGE BLANK NOT FILMED 
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APPENDIX A. MODEL FOR PLATFORM WITH FLEXIBLE CROSS ARMS 

This Appendix documents the application of standard finite element 
techniques to modeling the LSST reference platform with experiments rigidly 
attached to the platform arms. The platform arms are modeled as flexible 
appendages to a rigid bus. The solar panels are assumed to be rigid. 

A slx~degree-of»freedom model Is developed for planar motion of the platform 
cross arms. The equations of motion are developed using a lumped mass 
approximation and a consistent mass matrix approach. The double dlagonall- 
zatlon procedure for obtaining mode shapes and frequencies is discussed 
and the effect of parameter variations on mode shapes and frequencies Is 
illustrated. 

The purpose of this appendix is to provide a building block to more 
complete models of the LSST reference platform as developed In Appendices 
B and C. It is felt that this building block approach to modeling provides 
Increased understanding of the structural model which Is useful for controls 
design work. 

1. Configuration 

Figure A-1 shows the simplified configuration used for modeling of 
the platform cross arms. Motion is constrained to the pl 2 me. The motion 
of Interest for the initial modeling activity will be for rotations about 
the X axis. These rotations are tightly coupled with displacements along 
the z axis. For small angles the y axis motion is considered negligible. 
Hence, the model has six degrees of freedom, the vertical displacements 
^1’ ^^2’ ^^3 rotations 0j^, ^ 2 * The masses and (taken to be 

point masses) and Inertias I^ and I^ are associated with two rigidly 
attached science and applications packages. The mass M 2 and Inertia I 2 are 
associated with the central bus. The connecting elements are taken to be 
beams with length L, mass per unit length m, and flexural rigidity £1. 

2. Stiffness Matrix 

The stiffness matrix for a beam element is the standard one used for 
finite element techniques (see for example Ref. A-1, p. 158). With 
reference to Fig. A-2 we have. 


m m 


m 

6 -6 3L 3L 


1^ ■ 


. 2EI 

-6 6 -3L -3L 


^2 

i 






L 

3L -3L 2L^ 


®1 

T, 


3L -3L 2L^ 


X 


Here, and V2 are the vertical displacements of the beam end points, 
and and 02 are the rotations at these points (positive slope - positive 
rotation) . » ^2 corresponding applied forces and , T2 are 

the applied torques. 


Equation (A- 1 ) can also be written 


'fr 


6 3L -6 3l' 


\ 


« 2EI 

3L 2L^ -3L 


‘1 





I 

^2 

L 

-6 -3L 6 -3L 


'2 

^2 


Ji -3L 2L^^ 




Next, If we have two beam elements, we can find the combined stiffness 
matrix using tha standard direct stiffness approach. With reference to 
Fig. A -3 we have. 


^2 


6 3L -6 3l” 


h 

T 2 

. 2EI 

3L 2L^ -3L 


*2 





c 

^3 


-6 -3L 6 -3L 


*3 

1 

—1 

1 *** ‘ 


3L -3L 2L^ 


*3 




b. 






So that combining Eq. (A-2) and Eq. (A-3) we obtain 
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* m 

n 



^2 

. 2EI 



h 





re 

3L 

-6 

3L 

0 

0 

3L 

2L^ 

-3L 

L* 

0 

0 

-6 

-3L 

12 

d ’ 

-6 3L 

3L 

L* 

0 

4L* 

-3L 

L* 

0 

0 

-e” 

-3^* 

6 

-3L 

0 

0 

3L 

L* 

-3L 2L^ 


(A-4) 


Equation (A-4) defines the stiffness matrix (K) for the system of Fig. A-1 


3. Mass Matrix 

The consistent mass matrix is the standard used for finite element 
techniques (see Ref. A-1, p. 163 for example). 

If the mass per unit length of the beam in Fig, A-2 is m then the 
consistent mass matrix for this beam element is 


1 


ml 

WS 



’l56 

54 

22L 

-13L 



54 

156 

13L 

-22L 



22L 

13L 

4l2 

-3L^ 



-13L 

-22L 

-3L^ 

4l2 

» 

can be rewritten 





'156 

22L 

54 - 

13L 

ml 


22L 

4l2 

13L - 

3L^ 


54 

-13L 

156 - 

22L 



-13L 

-31^ 

-22L 

4l2 


'1 


(A-5) 


(A-6) 


The consistent Bass niatrlx for the second Bass eleBsnt is identical 
(since we assume m and L ere the same) , so the combined result Is 


156 22L 54 -13L 

22L 4L^ 13L -3L^ 

54 13L 312 0 


0 V 


0 0 e 

54 -13L V 


To the mass matrix of Eq. (A-7) we must add the discrete masses and 
Inertias of Fig. A-1 


so that the mass matrix for the system of Fig. A-1 Is 

with Mg taken from Eq . (A-7) . 

4. The Equation of Motion 


The equation of motion for the system of Fig. A-1 is 
M V + K V » F 

where 

T 

V - V202 V303I ■ nodal coordinate vector 

T 

F » tFj^Tj^ F2T2 F3T3] • force vector 

and M and K are as defined above. 





43 


In many cases, the effect of is negligible and the equation of 
motion may be approximated as 


Mjj V + K V - F 

with Mq defined by Eq . (A-8 ) . 


(A-10) 


5 , The Eigenvalue Problem 


If we set V ■ V f where v is a scalar and f is a vector of 
dimension 6, then Eq. (A-9) gives us 

(-8^ M f + K f) V e^®^ » F 
2 

now set F ■ 0 and s « X, so 


(A-11) 


X M f - K f 


(A-12) 


To put this in the standard eigenvalue problem form we will use a double 
dlagonalization procedure; 

First we chose 4^ so that 


F = g and M ® 1 


(A-13) 


here is the eigenvector matrix of M, with the individual eigenvectors 
scaled to satisfy Eq. (A-13). 

Applying Eq. (A-13) to Eq. (A-12) we obtain 


X g - K g 


(A-14) 


A g - K g 


(A-15) 


This is now in standard eigenvalue form. Now set 


g ■ *2*' 

where is the eigenvalue matrix for K. Then 


X h - *2 *2 ^ 


(A-16) 
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* 111 A 

Now, R la ayanetrlc and hence K la too. So ^2 K ^2 la diagonal 
and conalata of the elgenvaluea of the eyatem. The elgenvectora are 

hj^ - [1 0 0 0 0 0] 

h2 - [0 1 0 0 0 0] (A-I7) 

e 

h^ - [0 0 0 0 0 1] 

Theae can be tranaformed back Into the coordlnatea of Fig. (A-2) with 

♦l “ *1 *2 ‘‘i ^ 

This aolutlon procedure for the eigenvalue problem la well known and avoids 
Inverting the mass matrix. It also allows the eigenvalue routines to work 
with symmetric matrices throughout. Notice, however, that If Eq. (A>10) 

Is to be used, Is already diagonal and hence easily Inverted. In this 
case, there Is no need to use the double dlagonallzatlon procedure. 

The of Eq. (A-18) are the system mode shapes and the natural 
frequencies are given by 

“i * 1 - 1,...,6 (A-19) 

6. Efr.ccts of Parameter Variations 

Table A-1 presents eight sets of parameter values which were used to 
represent the platform arm model for Fig. A-1. These values are believed 
to be representative of a power system/platform similar to the 25 -kW 
power system reference concept proposed by Marshall Space Flight Center. 

Table A-2 presents the modal frequencies associated with the parameter 
values of Table A-1. The mode shapes are shown in Fig. A-& (the mode 
shapes do not vary qualitatively with changes In parameter values) . The 
following observations can be made;. 

• A comparison of the results for data sets 7 and 8 show that the 
mass per unit length of the beam element has a negligible effect 
on the frequencies. Hence, the consistent mass matrix approach may 
be set aside In favor of the simpler diagonal lumped mass matrix. 
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• A comparison of results for data sets 2 and 3 shows that the inertias 

(Il» l2* ^3^ have little effect on the first synnetrlc bending 
mode (mode 3) . 

• A comparison of results for data seta 4 and 5 shows that the inertia 

I2 has little effect on the second bending modes (modes 5 and 6) . 

• Increasing L, or reduces the modal frequencies and increasing 

El Increases them. 


^ . Comparison of Results to those for a Simple End Loaded Cantilever 
Beam 

The modal frequencies for the first bending mode may be compared to 
those of a simple end loaded cantilever beam (see Fig. A-5) . 

The tip deflection is (see for example. Ref. A-2, p. 518) 

^ " 3EI 

So the stiffness is 
V F 3EI 


The differential equation is: 

M y + K y «* 0 
y + I y - 0 

So that the natural frequency of vibration is Just 



K 

m 


3EI 


(A-20) 


using the parameter values of data set 2 we see 




3 X 38.45 X 10^ 
(3x10^) X (3)^ 


37.74 rad/sec 


6.0 Hz 
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Thla is cloaa to th« results of date sets 1-5 (5.3 to 7.3 Hz). Hence, to 
a first epproxlfflstlon, Bq. (A-20) can be used to estlaste the lowest 
frequency of the aysten of Fig. A-1. 


Table A-1. Pareneter Values 


Data i L 
Set I (a) 


El 2 
(Nhd^) 

xlO^) 

(*<8/m) 

(kg- 

xl03) 

M2 

(kg- 

xl03) 

M3 

(kg«2 

xl03) 

(kgn^ 

xl03) 

l2 

(kgmZ 

xlOh 

^3 

(kgn^ 

xlO^) 

38.45 

12.3 

3 

13 

3 

0 

86 

0 


0 

t 

i 

1 

i 

0 

0 

0 




Table A-2. Modal Frequencies (Hz) 












A10DE 1: RIGID BODY TRANSLATION 


MODE 2: RIGID BODY ROTATION 


MODE 3; 1st SYMMETRIC BENDING 


MODE 4: 1st ANTISYMMETRIC BENDING 


MODE 5: 2nd SYMMETRIC BENDING 


MODE 6: 2nd ANTISYMMETRIC BENDING 


Fig, A-4. 
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APPENDIX B. MODEL FOR PLATFORM WITH RIGIDLY ATTACHED EXPERIMENTS 

This Appendix documents the application of standard techniques to 
modeling the LSST platform with experiments attached rigidly to the 
platform arms. The platform arms and solar panels are modeled as flexible 
appendages to a rigid bus. 

Nine-degree- and eleven -degree-of -freedom models are developed by 
adding flexible solar panels to the model developed in Appendix A. 

1 . Configuration 


Figure B-1 shows a simplified configuration used for developing an 
eleven-degree -of -freedom model for the LSST reference platform. The two 
degrees of freedom 6^ and 6^ are later eliminated (see sections which 
follow) to obtain a nine -degree -of -freedom model. The variables 6,, 6<;, 0, , 
Vj, and are associated with the solar panels. The central bus mass 
and inertia are represented by M2, I2, and I5 (M2 * = 1/2 total bus 

mass, I2 * I5 = 1/2 total bus inertia). The two experiment (payload) 
packages are represented by M^, Ij^ and I3. The masses and M3 are 
assumed to be connected to the bus by massless beams of length and 
flexural rigidity EI^. The solar panels are represented by M^, and 
^6* ^6 ^^A ~ ^6’ ^A "" ^6^' masses M^ and M^ are assumed to be 

connected to the central bus by massless beams of length Lj^ and flexural 
rigidity EI^^. Since the bus is rigid we have 62 = 63. 


2. Solar °anel Model 


The solar panel model has six degrees of freedom (before applying the 
constraint that 62 = 65), The stiffness matrix found by finite element 
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A, 

Eq. 1 

(A-A)) 

is 

-6 

3L 

0 

0 

-3L 

l2 

0 

0 

12 

0 

-6 

3L 

N 


-3L 




. 6 

-3L 


2L‘ 



'' 4 ' 


e. 


4 




“5 



u 

/ 

. M 


symmetric 


Here F , i » 4,5,6 are the forces associated with the displacements 
i - 4,5,6; and 1 - 4.5,6 are the torques associated with the rotations 

e^, i - 4,5,6. 

A diagonal mass matrix is used (see Appendix A, Eq. (A-8)). 




I 


I 



3 . Constraint Impleme ntation (9^ = 

The stiffness end mass matrices for the platform have an identical form 
(Appendix A Eqs. (A-4) and (A-8». In order to combine them »e set - 15. 
and e, = 65 to obtain an 11 by llmassmattix and stiffness matrix. 

This procedure Is straightforward. First reorder the rows and columns 
of the matrices so that tor the platform arms we have the 6^ equation 


on the bottom: 



% 


' 

reordered 



stiffness 



matrix 


''2 

(or reordered 


''3 

mass matrix) 





®2 

W m 


For the solar panels reorder the tows and columns of the matrices so that 
the Tj, e, equation is on the top (and then set I 5 ■ ®5 “ ' 



— • 

reordered 


m 

®2 

stiffness 


V. 

matrix 


*4 

(or reordered 


®4 

mass matrix) 


'^5 

lb « 

. 

*6 






lU'IV'lUWJP 
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Now, combine these equations to obtain an 11 x 11 matrix. For 


example, the mass matrix becomes 


^ -1 

r 1 




"1 




h 


f2 


H 2 




•<3 


T3 


*3 


^2 

= 

‘'2"'5> 


f4 


N 4 


T 4 


>4 


fs 


"5 


fe 


«6 




•e 
' ^ 

1 

i. 


The stiffness matrix is similar, it will look something like 


Reordered . 

K matrix 

for platform ' 

arms r 

= 'cj 

" ” T Reordered 

I K matrix 
1 for solar 

I panels [. . 

► J ^ ^ 

The common element (C) will be the sum of two terms, one from each of the 
two stiffness matrices. 


4. Static Condensation 

What we have now is an eleven degree -of -freedom model. However, the 
degrees of freedom associated with the rotation angles 6^ and 0^ are not 
essential for a preliminary model. Excluding these degrees of freedom 
still provides two flexible modes for the solar panel (first symmetric 
and antisymmetric bending modes) as will be shown in what follows. 
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The process of eliminating 0^ and 0^ from the lumped mass matrix Is^ 
simply to set “ 0 and compress the matrix to a diagonal matrix (M) 
of dimension 9. 


0 and 0, can be eliminated from the stiffness matrix using a process 
4 o 

known as static condensation (see, for example. Ref, A-1, p. 172). To 
accomplish this we begin by reordering the elements of the stiffness matrix 
and partitioning it so that it has the following form 



or 


m m 



1 


• e 



'‘tt 

CD 






1 

<D 


1 

<x> 

> 

1 



F, - 

Vt . 


IVel' 


where 
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and *^ 0 t’ ^ 60 * P®*^tltioned elements of the reordered 

stiffness matrix. 


Now set ■ Ig ■ 0, Also set ■ 0, since we assume that no 

external torques will be applied to the masses 


*• « 



- 





K ^ 

K 


V 

t 

m 

tt 

t 0 



0 





V 



0 t 

00 _ 


0 _ 


then 


“et 't ■'se'e ' “ 


or 


\ - - ■'as'' ■'at ''t 


and so 


(B-1) 


^ ‘ ‘"tt - ''ta ■'ea'^ •‘at) \ 


So the reduced stiffness matrix is represented by: 


K = 


tt 


K 


te 




(B-2) 


Using K and M we can now work the problem with nine degrees of freedom 
instead of eleven, Also, Eq. (B-1) can be used to solve for 0^ and 6 ^ 
under the assumption that I 4 “ = 0. Notice that 0^ and 0^ will not 
be zero. 


5. The Equation of Motion 

The equation of motion for the system of Fig. B-1 is 
M Vt + K Vt * Ft 

A A 

where M, K, V^, and F^ are as previously defined. 

6 . The Eigenvalue Problem 

The eigenvalue problem has the same form as that given in Appendix A, 
Eq. (A-12) 


A M f » K f 






~i T"* *•' 
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2 

where X Is the scalar eigenvalue (^ * **>^ * natural frequency squared) and 
f Is the eigenvector (or mode shape) associated with X, 


7 . Mode Shapes for System Elements j 

r 

In order to visualize the system mode shapes, it Is instructive to 
first examine the mode shapes of the Individual elements. Our model has ; 

two types of elements as shown In Fig. B-2. Type 1 has 4 degrees of freedom | 

giving two rigid body modes and two elastic modes. Type 2 has 3 degrees of ; 

freedom, giving two rigid body modes and one elastic mode. 

Figure B-3 shows the mode shapes for the type 1 element. The type 2 i 

element has the same first three mode shapes shown In Fig. B-3 but lacks 

the 4th mode. 

1 

I 

The first two mode shapes Involve no strain energy (since the beam i 

element connecting the end masses Is undeformed) and have u 0. The 4th | 

n i 

mode shape has a higher strain energy than the 3rd mode and also a higher 
natural frequency. 


8 . Sample Computer Run Results 

Four computer runs were made with the eleven-degree -of -freedom model and 
two runs with the nlne-degree-of -freedom model. Table B-1 lists the parameter 
values used and Table B-2 the resulting natural frequencies. 

Comparison of runs 4 and 5 show that the results are identical for 
the first nine modes. In other words, run nine confirms the fact that 
eliminating 6^ and 6^ Is equivalent to setting I^ and Ig to zero. 

Comparison of runs 2 and 3 shows that Increasing the mass and 
Inertia of the experiments Mj»l3) results in lower frequencies for 

the bending modes of the platform arms, but has little or no effect on the 
solar panel modes. 

Comparison of runs 1 and 2 shows that Increasing L and decreasing El 
lowers the natural frequencies as we would expect. 


As poi nted out in Appendix A, the lowest frequency Is roughly propor- 
tional to J EI/ML^ or 




/ 


El 

ML' 


1 


; C s constant 
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Between runs 1 end 2 this ratio for the platform arms changed by an amount 

- 10.0 

“n2 ^ 

This agrees well with the coo^uter run results for the 1st syometric mode: 


( 1 ) ^ 

— - 3. 91/. 546 


7.2 for platform arms 


“n2 


.356/. 035 ■ 10 for solar panels. 


Comparison of runs 2 and 4 shows what the result of neglecting 
and is on the modes retained. The platform arm modes are unaffected 
and the solar panel mode frequencies increase somewhat (compare modes 8 and 
9 for runs 2 and 4). When using the nine -degree -of -freedom model, the 
flexural rigidity (El) can be reduced somewhat as compared with the eleven- 
degree -of -freedom model if the same modal frequencies are desired for the 
modes 8 and 9. 


The mode shapes for the system of Fig. B-1 are sketched in Figs. B-4 
through B-6. The mode shapes do not differ qualitatively from run to run. 
The shapes are various combinations of the element shapes sketched in 
Fig. B-3. Table B-3 lists the mode shapes quantitatively for run 6. 


9. Conments on Interaxis Coupling 

The mode shapes of Figs. B-4 through B-6 demonstrate clearly that x 
and y axis rotations are tightly coupled (see Fig. B-1 for axis definition). 

For example, any input which excites the symmetric bending modes of the 

* 

platform arm will cause the rotations 6^ and 6^ . But, because the symmetric 
bending mode results in a displacement V 2 , this motion will cause rotations 
of M 2 and about their y axes. This in turn will excite torsional 
modes (not modeled in this paper) associated with the solar panels and 
platform arms. In this way, we can see that all six masses will have 
rotations about their y axes as a result of an excitation of a platform 
arm symmetric mode. 


* 


The torques T 


1 


and T^t for example, will excite these modes (but T 2 will not). 
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The magnitude of this coupling will depend on the parameters of the 
model and may or may not be small enough to permit meaningful single axis 
control system designs. This coupling effect will be examined more 
carefully in future work. 
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Table B-i Run Parameters 


Parameter 

Run 1 

Run 2 

Run 3 

Run 4 

Run S 

Run 6 

Platform Arms 







La (m) 

5 

15 

— 



15 

EIa(NVxlO®) 

40 

20 

— 



20 

M, (kg X 10^) 

3 

3 

12 

3 

3 

3 

I, (kg-m^xlO^) 

1 

1 

4 

1 

1 

1 

M2 (kg X 10^) 

6 

— 




6 

I2 (kg-AlO^) 

60 

— 




50 

M3 (kg X 10^) 

3 

3 

13 

3 

3 

3 

I3 (kg.n?xlO^) 

1 

1 

4 

1 

1 

1 

Solar Panels 







Lj, (m) 

20 

— 




20 

EIj, (nVxIO®) 

10 

.1 




.1 

M4 (kg X 10^) 

0.5 

— 




.5 

I4 (kg-m^xlO^) 

50 

50 

50 

0.001 

0 

0 

Mj (kg X 10^) 

6 

— 




6 


X 

X 

X 

X 

X 

X 

Me (kg X 10^) 

0.5 

— 




.5 

Ig (kg.m^xlO^) 

SO 

SO 

SO 

0.001 

0 

0 



Table B-2 Natural Frequencies (Hz) 


Rigid Body Modes 

Run 1 

1 (solar panel translation) 

0 

2 (arm translation) 

0 

3 (rotation) 

0 

Platform Arm Modes 

4 (1st symmetric) 

3.91 

5 (1st antlsynsnetric) 

5.75 

6 (2nd symmetric) 

29.3 

7 (2nd antisymmetric) 

29.3 

Solar Panel Modes 

8 (1st symmetric) 

.356 

9 (1st antisynmetrlc) 

.584 

10 (2nd symmetric) 

1.33 

11 (2nd antisymmetric) 

1.37 


Run Z Run 3 Run 4 Run 5 Run 6 

» 0 

0 

0 


.546 

,'•30 

.546 

.546 

.546 

2.04 

1.98 

2.04 

2.04 

2.04 

11.7 

5.86 

11.7 

11.7 

11.7 

11.7 

5.95 

11.7 

11.7 

11.7 


.035 

.035 

.047 

.047 

.0471 

.039 

.035 

.049 

.049 

.0493 

.133 

.133 

22.5 

X 

X 

.130 

.129 

22.5 

X 

X 
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Table B-3 Mode Shapes for Run 6 


Mode Number and Displacement 


1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

-9JE-3 

-l.lE-2 

9.1E»3 

2.4E-3 

7.4E-4 

7. BE -4 

0 

-6.0E-3 

0 

-9.1E-3 

0 

-9.1E.3 

0 

-7.4E-4 

0 

0 

0 

0 

-9.1E-3 

l.lE-2 

9.1E-3 

• 

-2.4E-3 

7.4E-4 

-7. BE -4 

0 

6.0E-3 

0 

0 

7.4E-4 

-USE -3 

.2.5E-3 

2.2E-2 

2.2E-2 

0 

4.0E.4 

0 

0 

7.4E-4 

1 .8E-3 

-2.5E-3 

-2.2E2 

2.2E-2 

0 

4.0E-4 

0 

0 

7.4E-4 

0 

4. IE- 3 

0 

4.6E-4 

0 

3.9E-4 

1 ,2£-2 

0 

-1.5E-2 

0 

4.0E-5 

0 

1 .3E-7 

2.9E-2 

•2.8E-2 

1.2E-2 

0 

0 

0 

0 

0 

0 

-4.9E-3 

0 

1 .2E-2 

0 

1.5E-2 

0 

-4.0E-5 

0 

-1.3E-7 

2.9E-2 

-2.8E-2 











MODE 1: RIGID BODY TRANSLATION 


MODE 2: 


A10DE 3: 


MODE 4; 


RIGID BODY ROTATION 


1st BENDING MODE 


2nd BENDING MODE 


C 


C 


C 


B-3. Mode Shapes for Type 1 Element 



MODE 1: 

RIGID BODY TRANSLATION OF 
SOLAR PANELS 

o 

MODE 2: 

RIGID BODY TRANSLATION OF 

o 


PLATFORM ARMS 

o 

MODE 3: 

RIGID BODY ROTATION 

6 



o 


Fig. B-4. 

Rigid Body Modes 



MODE 4: 1st SYMMETRIC BENDING MODE 



MODE 5: 1st ANTISYMMETRIC BENDING MODE 



MODE 6: 2nd SYMMETRIC BENDING MODE 


MODE 7; 2nd ANTISYMMETRIC BENDING MODE 





Fig. B-5. 


Elastic Modes for Platform Arms 
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MODE 8: 1st SYMMETRIC BENDING MODE 



a 




MODE 9: 1st ANTISYMMETRIC BENDING MODE 



MODE 10: 2nd SYMMETRIC BENDING MODE 



MODE 11: 2nd ANTiSYMMORIC BENDING MODE 



Fig. B-6 


Elastic Modes for Solar Panels 


APPENDIX C. TRANSFER FUNCTIONS 

This Appendix documents transfer functions for the nine -degree -of - 
freedom model of Appendix B. Transfer functions are developed from modal 
data. Comments on observability and controllability are made for torque 
actuators and angular position sensors on the central bus or on the 
platform arms. Transfer functions are included for noncolocated actuators 
and sensors. 

1, Transfer Function Derivation 

This section addresses the derivation of transfer functions between 
a force (or torque) at any point in a structure to a sensor at any point 
in the structure. 

For this purpose the system model is assumed to have n degrees of 
freedom and be of the form 

Mx + Kx = F = Bu 
y = Cx 

where 

X is the n "dimensional" nodal coordinate vector 
u is them "dimensional" control input vector 
F is the n "dimensional" force vector 
y is the H "dimensional" output vector 
M is the n X n mass matrix 
K is the n X n stiffness matrix 
B is the n X m control distribution matrix 
C is the £ X n output matrix 

Let t be the n x n matrix of eigenvectors for this system such that: 

M = I 

and 

T 2 

4> K ♦ = A 

2 

where I is the identity matrix and A is the diagonal matrix of eigenvalues 
(or modal frequencies squared) . 
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Then Eq. (C-1) can be rewritten In terms of the modal coordinates q 
(x ■ ♦q) as 

“2 T 

q + A q - ♦ F (C_ 2 ) 

Suppose we now assume modal damping* and write 

** • A *T* 

q+Dq+A‘q»# F (C-3) 

where D is a diagonal damping matrix. 


Now take the Laplace transform of Eq. (C-3) to obtain: 
q(s) = (s^I + sD + A^)"^ F 
or in terms of the original coordinates 
x(s) = <Ks^I + sD + A^)"^ F(s) 


Equation (C-5) gives the transfer matrix relating F to x . The transfer 
matrix between u and y is given by 

y(s) = C *(s^I + sD + A^)~^ 4 ’’ B u(s) (C-6 

To find the transfer function between any force (^j) ^od nodal coordinate 
(Xj^) we proceed as follows: define 

Q(s) = (s^I + sD + A^)“^ 

= diag. (Q^(s) Q 2 (s) ... Qj^(s)) 

Q^(s) = (s^ + 2 u)^ + 

<<>2 . • . 4 ' jj ] 

and 

*^i “ ^“^li '*’2i ••• *^ni^^ 


A 

See, for example. Ref. A-1, pages 194-199 for a discussion of this assumption. 
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where 

Is the Jth component of the 1th eigenvector 
ith modal damping coefficient 
■ ith natural frequency 

Then, from Eq. (C-5) we have 


X(s) - (4i^ (>2 ... 

Ql(s) 


Fa ’ 
*1 

1 

F(s) 

1 


Q2(s) 

• 



1 

! 

1 


• 

■ Q„(s) 


A T 
_^n 


n 

= ( I Q.(s)4*/) 

i-1 ^ 1 1 

F(s) 


(C-7) 


Equation (C-7) is an alternate representation of the transfer matrix given 
in Eq. (C-5). Now define 


x(s) = 

/ ^T 

'^1 ^2 * * * ^n. 

F(s) = 

(Fj F2 ... 

that 



n 

x^(s) • 

i»l 


(C-8) 


Equation (C-8) defines the transfer function between the force F. and the 

J 

nodal coordinate x, . 

k 


Using (C-7) we can write (C-6) as 
n 

y(s) - C( I 4>. Q. ) B u(s) 
i«l 1 1 l 


( ^ C ♦ Q ^ B) u(s) 
i-1 


(C-9) 


now set 


B - [B^^ B2'‘^ ... bJi"^ 


C “ [C. Qj C£,] 
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where each Is a 1 x m row vector and each is a n x 1 column vector. 
Then 


and 


“ 1 X m row vector 

c 4 

» £ X 1 column vector 
n n 


y(s) » [ I ( I C ^ ) Q (s)( I <?, B )] u(s) 

i»l j-1 J ^ j=l J 


i I C Q (s) B,] u(s) 
i«l ^ ^ 1 


(C-10) 

(C-11) 


now, define 


®i ' ^®il ®i2 ••• 


" ^^li Si ••• Si^ 


where and are scalars. Define 

y(s) = (yj^ y 2 ... y^^)^ 

u(s) = (u u ... u 
X / m 

so that ^ 

y^(s) = Qi(s) B^j] Uj(s) (C-12) 


Equation (C-12) defines the transfer function between the control input 
and the output 


u 


j 


2. Transfer Function Symmetry 

Equation (C— 8) demonstrates that the following two transfer functions 
are equivalent: 
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This r..uU 1. known a. Maxwell., i„ „f d.fi.„i<„, 

In 1864) . Ih, ajoBatt, aaan hare arlaaa bacauaa ol the aynnatry of the 
original M and K matrices , 


3. Transfer Function Zeros 

The nurt.ar of «ro. (Nz) In .h. nuoorator of Eq. (c-7) (or (c-12)) 1. not 
at all obvlona. An uppar bound la »z . n . 2. Thl. can bo soon by expanding 
tq. (C-7) over a common denominator which Is the product 

n 2 

X (s^ + 2 ; u). s + 

1-1 ^ ^ i 

However. In oany oaaes the degree of the nueerator polynonaal will be le.e 
than n-2 when k?<J . 


As an example, consider the system of Fl«. C-1. For this system the 
transfer function Xj(s)/F^(s) has 6 zeros. X 2 vs)/F.(s) has 4 zeros, 
x3(s)/F^(s) has 2 zeros, and x^(8)/F^(s) has no zeros. 


This can be demonstrated as follows 

(Kij ^ + Kj) 


The system equations are 



To find the inverse of the matriv 

tne matrix (Ms +K) we can use Cramer's rule 

numerator of the Inverse has the form 


The 


* — — 

See for example. Ref. C-1, p. 494. 
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■u 

"12 

"l 3 

■« 


■22 

■23 

■m 

■m 

*32 

■33 

*>4 

*« 

"42 

"43 

■4* 


and the denominator is the determinant of (Ms^ + K). Each element la 
found by deleting the ith row and Jth column of Ms^ + k) and the 

determinant with proper sign. So 


N 


'll 




’"2 

0 


*^3 ”*3 




In this case, the coefficient of the s® term ia M2 so N.^ has degree 


6 . For Nj^ 2 we have: 


Kl -«2 


"12 


“ 3«^3 -«3 


0 




In this case, there is no s® term and the s^ term coefficient is K, M, M, . 

1 3 4 


For we have: 


’ 13 ’ I ® -*2 

0 0 






. E: ■ ■■■ ■ II ■ - 
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So, in chla cmo, %m havo o^ aa cha hlghaat dagraa of a. 



For N 


wa have: 



-«2 

0 -Kj 


0 0 





80 ia juat a constant tarn, and there are no zeroa in the tranafer 
functiona x^(8)/F^(a) - Xj(8)/F^(s). 

So we conclude that for chain systems of this type we can determine 
pretty much by inspection the number of (finite) zeros. For other systems 
this is not so easy, especially when the M and K matrices are full or nearly 
full. 

It should be noted, however, that if zeros at infinity are Included, 
then every transfer function can be thought of as having the same number of 
zeros as poles. 

Suppose that we compute the transfer functions using Eq. (C-7) . Then 
because of computational inaccuracies the order of the numerator polynomial 
for each transfer function will almost always turn out to be n-2 for systems 
of any size. 

If a polynomial root finder is employed it will attempt to find n-2 
roots. It may or may not be successful and if it is successful, some 
of the roots found may have very large magnitudes. The user of such a 
program should understand that the most probable explanation for the root 
finder bombing out or obtaining roots with very large magnitude is that the 
true system has fewer zeros than n-2. 

If the root finder does bomb out, it is necessary to reduce the 
polynomial to a lower degree (by sisqile truncation) before attempting 
another solution. If the root finder obtains very large values for the 
zeros (as compared to the pole magnitudes) then these zeros should usually 
be discarded (which is equivalent to replacing them by zeros at infinity). 
Some Judgement is required, however, on what is “too large." 
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4 , for Nln«" 0 «Rr«t-of -^ retdoa Mod«l 


Th« initial control problM for thU structure is related to the angles, 
02 , 03 where 0 ^ and 0 j are the experlaent pointing angles and 02 i* the 
central bus orientation angle. 


The associated control torques are defined as T^, T 2 , and T 3 ; where 
T^, 1 • 1 , 2,3 is the torque applied at I - 1,2,3 about the x axis. 

The transfer functions of Interest are: 


( 8 ) 


transfer function between a torque applied at and the 
rotation angle at same location. 


®2 

^ ( 8 ) 

2 


transfer function between a torque applied at the central 
bus and the rotation angle at same location. 


^ ( 8 ) 


transfer function between a torque applied at the central 
bus «id the rotation angle at 


(s) 


transfer function between a torque applied at and the 
rotation angle at 


Also of interest are the following transfer functions 


and 


®1 

(8) - (s) 

Ti I2 


®2 ®1 

=4 (8) - si (S) 

^1 ^3 


There are many other possible transfer functions of lesser interest 
which will not be discussed. 


Figures C-2 through C-5 present pole/zero configurations for the four 
principal transfer functions resulting from the model of Appendix B using 


the Run 6 date set parameters. For simplicity of platting no rfamping has 
been added. The gain (K) of the transfer function has also been shown using 
the following convention 

K (8-2^)(s-Z 2) ... (s-Z^) 

" "(s-Pji^s-Pj) ... (s-pj 

where 

= the 1th zero 
and 

Pj = the jth pole. 

5. Observability and Controllability 

The pole/zero plots of Figs. C-2 through C-5 give us information on 
observability and controllability. Whenever pole/zero cancellation occurs 
this is an indication that either observability or controllability (or both) 
has been lost (see next section) . 

e. 

Figure C-2 shows us that for ^ (s ) there is a pole/zero cancellation at 
w = 0.296 rad/s. This frequency corresponds to the first symmetric mode of 
the solar panels. Any motion observed at as a result of solar panel motion 

must be transmitted through the bus. The symmetric modes of the solar 
panels, however, are not transmitted to the platform arms (see Appendix B) 
so this indicates a lack of observability. Furthermore, this mode is not 
controllable since rotations and displacements of the platform arm cannot 
excite the symmetric modes of the solar panels. 

e. 

Figure C-3 shows us that for ^ (s) there are pole/zero cancellations at 
0 ) = 0,296, 3.43, and 73.3 rad/s. These frequencies correspond to the 
symmetric modes of the solar panels and platform. Since the symmetric modes 
do not affect 6^ these modes are not observable, and since a torque T^ 
cannot excite these modes they are not controllable either. 

Figure C-4 shows us that for (s) and ^ (s) there are pole zero 
cancellation at a> = 0.296, 3.43, a^ 73.3 rad/L As before the solar 
panel mode is neither controllable or observable. This is true for both 
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6 Oj 

(s) and s— (s). The platform arm symmetric modes are observable in (s) but 

^ ^1 ®2 ^2 

not controllable. In the case of ^ (s) the reverse is true: the 

platform arm symmetric modes are controllable but not observable. 

®1 ®3 

Figure C-S shows us that for (s) and ^ (s) all the modes are control- 

I3 i 

lable and observable except the solar panel symmetric modes which are 
neither controllable nor observable. 

It should be noted that controllability and observability of a particular 
mode are necessary if we are to change the dynamics of a system related to 
that mode, but do not in themselves assure us that a satisfactory control 
system design can be achieved. For example, the configurations of 
Figs. C-3 and C-4 will be much more difficult to control because of the 
consecutive pole patterns on (or near) the Imaginary axes, and because of the 
right -half plane zero(s). 

It should be noted that although the symmetric modes of the platform 
arms and solar panels are not controllable or observable at the central 
bus when considering 62 and T 2 , they become both controllable and observable 
when torques and motions about the y-axis of the bus are included. The 
present model does not include these torques explicitly , but they are 
implicit in the forces F 2 and F^ and the displacements V 2 and V^. 

6 . More on Pole Zero Cancellations 

In the previous section we used arguments associated with the mode shapes 
to determine whether pole zero cancellation arose from loss of controllability 
or from loss of observability, or both. In this section we will make the 
argument a bit more precise. 

The system equations as stated earlier are: 

Mx + Kx » F 
and 

i q • F (C-13) 

where 


X » ® q 


(C-14) 


and ^ is the eigenvector (mode shape) matrix. Here x represents physical 

coordinates, and q the modal coordinates. For illustrative purposes, 

T T 

consider a 3 -dimenslonal system with x ” [x^^ X£ x^] , q " [q^^ q2 qsl , 

F - (F^ F2 F3l'^, and 



*21 

♦31 

*22 

♦32 

*23 

*33 


T 

Here ^^2 ^1^3! represents the ith mode shape. Now, Eq. (C-13) can be 

written as follows 
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k m 
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“ 2 

^2 ^ "2 ^2 
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m 
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» m 


From this it is clear that for to be influenced by 
(this is Che controllability issue) . 


we must have ^ . 

mj 


* 0 


Similarly, Eq. (C-IA) can be written 
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*31 

*12 

*22 

*32 

*13 

*23 
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L 

^ m 

^1 

1 

‘*2 

) 

^*3 


From this we see that for q Co influence the output at x. we must have 

m K 

1* 0 (this is the observability issue). 


Now consider an n-dlmensional system. The transfer function between 
the Force F^ and output was found before (see Eq. (C-7)). For the 
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undasped case we have: 


f7<s) 


8 + 


2 . 

8 + U, 


+ . . , + 


2 ^ 2 
m 


+ . . . + 


2 . 2 
8 + W_ 

Xk 


(C-15) 


Now, if either ■ 0 or « 0, then the term involving 8 + 

disappears the transfer function has no pole at s^ + • 

Hence, if — (s) is written with the term s^ + u ^ in its denominator, it 
r j m 

must have the same term in the numerator. In other words, we must have pole 
zero cancellation. Hence, we can examine the entries in the eigenvector 
matrix to determine controllability and observability; or we can look for 
pole zero cancellations in the transfer functions and if they occur look 
at the eigenvector matrix to determine their cause. 

The discussion in this section and the last has implicitly assumed 
that the natural frequencies are distinct. If they are not distinct (for 
example we may have several rigid body modes) then our comments must be 
modified slightly. 


Looking at Eq. (C-15) we see that pole zero cancellation will certainly 
occur if “ 0 or « 0. However, if u>j^ * (for example) then it is 
possible to have pole zero cancellation even when and ^ 2 j ^' 21 ^ non- 
zero. Thus 1 0, (fr , 0 is necessary for controllability/observability 

mj mk 

but not sufficient. If, however, ♦ . ♦ , 0 and u is a distinct (nonrepeated) 

mj mK m 

eigenvalue, then this mode will be controllable and observable (♦^j 0 is 

necessary and sufficient). For additional discussion on this topic the reader 
is referred to Ref. 2-1. 
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APPENDIX D. MODEL FOR PLATFORM WITH TWO HINGED EXPERIMENTS ATTACHED 

This Appendix extends the model of Appendix B to include two hinged 
experiment psckages, each having a single rotational degree of freedom. It 
is assumed that these packages are rigid. The resulting model has eleven 
degrees of freedom. The mass and stiffness matrix for the combined system 
is given. 

1. Configuration 

Figure 2-1 shows the configuration for the Reference Platform. This 
configuration is the same as that in Appendix B except for the two experiment 
packages represented by the masses additional degrees of freedom 

are given by the angles and y^. Yj^ end y^ nre taken to positive for a 
rotation about the plus x axis. Yj^ is measured from the negative y axis and 
Y^ Is measured from the positive y axis (l.e. for the nominal configuration 
shown in Fig. 2-1. both Yj^ and Y 3 sre zero). 

2. Development of Equations of Motion for an Experiment Package 

Figure D-1 shows a sketch of the experiment packages located on the 
y axis of the platform arm. Figure D-2 shows a free body diagram. The 
model developed will be for Yj^ snd y^ near zero so that forces along the 
y axis will not be considered. Symbols are defined as follows: 

* torque applied about the x axis of the experiment package located 
on the plus y axis. 

F^^ ■ reaction force applied by the base mount to the hinge point (P^) 
of the experiment package in the direction of the z axis. 

> point of force application 

Y^ - angle of rotation about the x axis relative to the Inertially 
fixed axis y. 

* inertial displacement of point P^ along z axis. 

“ mass of experiment package at location 3. 

1 ^^ > inertia of experiment package about hinge point P^. 

« distance between hinge point P^ and center of mass of experiment. 

Ti, F^^, P^, Yj^, are defined similarly. 
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Nov, If the point P^ were Initially fixed, then we would have 
“^3 " ^e3”'^3 

But, since Is not fixed we must account for the base acceleration. The 
proper relationship in this case is (Ref. D-1, p. 146) 

M*p xMr*H (D-1) 

-P “C -p -P 

where 

- torque about point P^ 

- vector from point P^ to cm 

M “ mass ■ M , 
e3 

r^ > acceleration of point P^ relative to Inertial space 

rate of change of angular momentum relative to point P^. 

The model being considered will consider only motion along the z axis and 
rotation about the x axis (actually there is also rotation about the y 
axis because of x-y axis coupling but this effect will be neglected In 
what follows). 


For Yj near zero we have 


M - T-e 
-p 3 X 

P “L-(e + •y- e ) 

-c „e3 y ’3 z' 

r - VoC 

— p 3 Z 

H = I ® 

-p e3'3 X 


The assumption here Is that x is a principal axis for the experiment package 
and that base motion is along z on] 
the X, y, and z axes respectively. 


and that base motion is along z only e^, e^, and e^ are unit vectors along 


Substitution into Eq. (D-1) gives 


T, e„ “ L - (e + y-jC ) x M ,V_e » 1 ,y- e 
3 X e3 y 3 z e3 3 z e3M > 


(D-2) 


or 


^3 " ”e3^e3^3 * ^e3^3 


^^3 " ^3^e3'^3 ^e3^3 


(D-3) 
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Li 


Now, we can obtain the force equation by application of Newton Law. 


F ,e - M , X (acceleration of CM of experlnent package) 
v3 < e3 


Since the position of the CM of the experlnent package relative to Inertial 
apace la given by (Vj + 


'v3 ■ “.3 ‘^3 ♦ '■.3’'3> 




Equations (D-3) and (D-4) gave the force and torque applied to the 
experiment package at point Pj. The force and torque applied to the platform 
arm at this point Is the same with a minus sign. 


For the experiment pacl'age on the -y axis of the platform armEqs. 
(D-3) and (D-A) must be modified slightly. Equation (D-2) becomes 


Vx - ‘-el <-y - ".iVx ■ ‘elU ^ 


or 


- ’ Mel^^el^l + ^el^l 


(D-5) 

(D-6) 


Also, since the CM position for this experiment package is 
we have 


•"vl ■ «el - h'-el> 


<D-7) 


To summarize what we have so far: 

For the experiment package on +y axis (location 3) the forces and torque 
applied to the experiment package are 


^3 ■ «.3‘-e3'‘3 * ‘,3’'3 


"v3 ■ "e3‘^3 ♦ ‘-.ji'V 


For the experiment package on -y axis (location 1) the forces and torques 
applied to the experiment package are 


h * ‘.I’l 




The forces and torques applied to the bus are of opposite sign. 
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3* Tht Mm« Matrix 

Defina the new state variable vector to be 

(ViVjVj OjSjej VjVjV, TjYjl'' 

The equations of interest are: 


'1 • 


- fsl - “.1 <'l - '■.A* 

(D-8) 



■ ',3 ■ '*•3 ^^3 '■«3’'3* 

(D-9) 

■ ^1 

(D-10) 

fs • '' 5 . - ^3 

(D-11) 

^10 ■ - ".I'-el'l * ■ ■'l 

(D-12) 

"a - «.3'-.3^3 * hlh ■ ■'3 

(D-13) 


where 

F, - sum of external forces (F ,) and experiment package reaction 
1 al 

forces (F^^ ■ - F^j^) applied to 

F^ ■ sum of external forces (^^3) experiment package reaction 
forces (F^j ■ - F^^) applied to Mj. 

F^ ■ sum of external torques (F^^) and experiment package reaction 
torques (-Tj^) applied to Ij^. 

Fj - sum of external torques <F^j) and experiment package reaction 
torques (-T3) applied to 1^. 

Fj^q ■ ^3 “ torque applied to experiment package at location 3 about 
hinge point P^. 

Now, let m^j 1 ■ 1,9 
j - 1,9 

be the mass matrix elements of Appendix B. Then, by using Eqs. (D-8) to 
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(D-13) we can obtain the new oass natrlx ee 


“33 " ®33 “*■ **03 

*1.10 “-“.l^-ei 

“ 3,11 " “e3^«3 


*10,10 

* 


*11,11 

-'e3 


“io.i 

■ *1,10 


“n,3 

• “ 3,11 


®10,j 

• 0 

for 

“ii.j 

■ 0 

for 

“i.lG 

• 0 

for 

®l,ii 

- 0 

for 

“ij “ 

®ij 

for 


all J except 1 and 10 
all j except 3 and 11 

all i except 1 and 10 
all 1 except 3 and 11 

all other 1 and j . 


4. The Sttffneai Matrix 

The stiffness matrix is unchanged frmn that of Appendix B except for 
the addition of xero elements to Increase the dimension from 9 to 11. This 
is because the rotations and Y 3 caii be made without any strain energy 
(i.e. there are no forces or torques proportional to or y^). 
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APPENDIX B. A TBCHNIQtlB FOR RB8BAPIM6 THE RIGID BODY MODES OF A SPACECRAFT 

This Appendix presents s techniqua idiich csn be used to reshape the 
eigenvectors associated with a coanon eigenvalue. This is of particular 
Interest for application to reshaping rigid body node shapes for flexible 
spacecraft. Standard eigenvalue/eigenvector programs do not necessarily 
give rigid body mode shapes which are pleasing. The algorithm presented 
here allows considerable freedom for reshaping the modes. 


(E-1) 
(t:-2) 

- diag. (Xj, X], . . . 

1^ ■ Identity matrix 

♦ Is the matrix of system eigenvectors, or mode shapes. A la the matrix 
of system natural frequcn<.les. 


1. System Equations . 

The system to be considered Is of the form 

Mx -f Kx • 0 

where M Is the system mass matrix (symmetric , positive definite) 

K Is the stiffness matrix (symmetric, positive semidefinite) 
X Is the coordinate vector (dimension n) 

Suppose we set x ■ 4 q 
where 4 satisfies 
4^ M 4 - I 

n 

4^ K 4 « 


2. Reshaping Algorithm 

Suppose the first a eigenvalues are repeated (for example, all sero 
for the rigid body modes). 

s.t ♦ • [«1 . . . ♦„] 


• • • 


«lMr« U th« «lg«nv«etor. Th«n ati^poM wt with to vothtpt 
♦jt ♦ 2 * • • • ♦« • ptrtlcultr tuppoM • [^34 ^2J * * * 
and tui^tt vt with to tlinintto 4 ^ fot toat i tnd for otch J 

J - 2,3, 

In othtr vordt, we with only 4 ^ to contain the component i t 

4^1 ^0 , 4ij - 0 3 - 2,3 m . 

We can accomplith thia by proceeding at follows: 

First eliminate 4^^ from 4j • This can be accomplished by setting 





where a 


+ a 4 , 


"12 


'll 


now and 4j^ are no longer orthogonal with respect to M and K so we 

must modify 4^ . So set - 4^ + b 43 • Then to maintain the orthogonality 

conditions of Eqs. (E-1) and CE-2) we must have 


SSfV* ♦ 

4i M 4j ■ 0 
4^ K 4j • 0 

4i> - 0 

« 0 


likewise we must have 

( 4 J + b 42 > K ( 4 j + a 4 ^) • 0 
or b 42 K 4j + • ® 

but if 4j K ^2 " • ^1 ***^ 


to we must have 

( 4 J + b 4j) M * * 

or a 4 j M 4j^ + b 4 ^ ** ^2 


to that b 


4 . M 4i 

- - a-= -a 


M 4, 


(E-3) 

(E-4) 


Now note that 


+ b ♦J) M ( b ♦j) - M ♦^ + “ ^2 

• 1 + b^ - 1 + e^ 

end (dj + a dj ) M + a ♦^>-♦2’^ M ♦j + ♦x^ “ *1 

. 2 

• 1 + a 

So to ranomallze and ^ wa must divide each of them by-Jl + a^. 

This results in the following: 

♦x •<♦!“* + a^) 

♦2 “ (^2 + « + a^) 

^i2 

where a » - -7= 

♦il 

This procedure can now be repeated using and ^3 to eliminate the 
coiq>onent 41^3 from ^3 , and so on until we have a new set of m rigid 
body eigenvectors, only one of which has the component t ^ 0 . 

Now we can work with m -1 eigenvectors, and eliminate a different 
coordinate from all but one of these; and so on until for the m^** 
eigenvector m -1 components have been eliminated. This procedure is 
Illustrated in Fig. E -1 for the cese M • 5 « Note that at each step any 
one coordinate can be eliminated from the remaining eigenvectors. These 
coordinates need not be consecutive as shown in Fig. E- 1 . 

3 . A Property of the Reshaping Algorithm 

Suppose that is to be used to eliminate 41^2 from • Then we can 

A 

show that bas the same magnitude as . 




In partleular wt havn 



•, -'twy -tf 


•t 


This property shows ue that if begins with s Urge velue for 
then It will still hsve this Isrge velue efter the resheplng process 

In other words, we ere assured of a car tain aaount of coq>utatlonal 
stability. 


Examples of Ways In Which the Reshaping Algorithm Can Be Used 

As an illustration pf how the reshaping algorithm might be applied 
consider the following example. Suppose the rigid body modes for a space- 
craft Include rotation and translation of the Center of Mass (CM). An 
attitude control system does not control CM translation. If there are 3 CM 
translational coordinates for the spacecraft, then it would be desirable to 
eliminate CM translation from all but 3 of the rigid body modes. This can 
be accomplished using the reshaping algorithm of this appendix provided that 
the eigenvectors contain coordinates defining the CM location, if the 
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eigenvectors do not contain the CM coordinates , then they can be augnented 
so that they do by using the mass properties of the spacecraft. 

If the CM translations are not eliminated from the eigenvectors used 
to design an attitude control system, then there will be uncontrollable 
modes associated with the design problem. These uncontrollable modes can be 
eliminated by reshaping followed by truncation. 

As a second example consider the following control problem. Suppose 
we wish to design an attitude control system for a single spacecraft axis 
using classical design techniques. If the spacecraft model has 3 rigid 
body rotational degrees of freedom then two of these are uncontrollable 
when performing a single-axis design. If the actuator /sensor pair used 
to accomplish single-axis control is colocated, then we can eliminate the 
rotational coordinate with which they are associated from all but one rigid 
body mode using tlie reshaping algorithm. If they are not colocated we could 
choose either the coordinate associated with the actuator or the one associated 
with the sensor and eliminate it from all but one mode. In this way we can 
eliminate two rigid body rotational modes from the spacecraft model which 
are either uncontrollable, unobservable, or both uncontrollable and unobservable. 
Through this process we can achieve a cleaner mathematical system model for 
use in the control design problem. 
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Fig. E-1. Reshaping Process for Mode 
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APPENDIX F. TZMB RESPONSE PLOTS FOR RUM FIO 


In this Appendix time plot* am praaantad which llluatrata the tranelent 
reaponse for a repreaentatlve alnulatlon (Run FIO of Table 4-1). Slmulatlona 
wem conducted ualng ACSL (Advenced Contlnuoua Sleulatlon Language) on an 
1100/81 Unlvac computer. 

Figure F-1 ahowa the experiment 1 actuator torque and angular reaponae. 

Flgture F-2 ahowa the central bua actuator torque and angular reaponae. 

The actuator torque and angular reaponse of experiment 3 la ahovm In 
Fig. F-3. Commenta on Individual plota follow. 

1. Figure F-1 

The actuator torque (Tj^) and angular reaponse (Yj^) for experiment 1 are 
shown in Fig. F-1. The commanded torques (T^j^) for this run were +20 N-m 
for the first 14 seconds, -20 M-m for the next 14 seconds, and zero r.hereafter. 
The commanded angle consistent with this. The actual torque 

applied shows the strong influence of the structural elastic response fed 
back to the controller by the rate and position sensor. The two dominant 
frequencies eeen in the torque curve are at .30 rad/s (.048 Hz) and 2.8 rad/s 
(.45 Hz). Note that the time plot results of Fig. F-1 agree with the closed 
loop eigenvalues plotted in Fig. 4-6. Figure 4-6 along with Table 2-2 show 
that the .3-rad/s frequency is associated with the solar panel asymmetric 
mode (open loop .4-rad/a) and that the 2.8-rad/s frequency is associated with 
the arm symmetric mode (open loop 5.7 rad/s). 

2. Figure F-2 

The actuator torque (T 2 ) and angular response (S^) are shown in Fig. F-2. 
The commanded angle (0^2^ aero. The primary response is a rigid body 
rotation which occurs because the bus controller frequency is low (0.01 Hz) 
compared to the frequency of the disturbance input. The only structural 
vibration frequency evident in this plot is at .048 Hz and is associated with 
the solar panel asymoietric mode. 
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3. Figure F-3 

The actuator torque (T^) and angular response (Y 3 ) are shown in 
Fig. F-3. The commanded angle (y^j) was aero. The frequency content of 
Tj is very similar to that of Tj^ for this run. In fact, after about AO 
seconds, the two torques are almost identical. Both of these torques show 
the influence of the structural vibrations associated with the solar 
panel asynanetric mode (.048 He) and the arm symmetric mode (,A5 Hz). Note 
that the oscillation caused by the solar panel asymmetric mode is + 1 prad 
at the 100-second point of this plot. This oscillation all by Itself is well in 
excess of the .05-urad pointing requirement presented in Table 3-1. 
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Fig. F-3. Actuator Torque and Angular Response for Experiment 3, Run F-10. 
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